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Abstract 

It is shown that the combination of the structure functions Fi pre- 
dicted by the CKMT model at low and moderate values of Q 2 and 
MRS99 parton distribution functions at high Q 2 gives a good descrip- 
tion of the data over complete measured region of x and Q 2 . Using 
these structure functions the main characteristics of the muon-nucleus 
inelastic scattering are calculated. Nuclear effects and contributions 
of the neutral current and 7 — Z interference are taken into account. 
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1 Introduction 



Muon inelastic scattering off nuclei plays an important role in muon propa- 
gation through matter. In this process muon may loss a significant part of its 
energy and can be scattered at large angle. Therefore, muon-nucleus inelastic 
scattering is of interest for numerous applications related to muon transport 
in matter, in particular, for calculations of muon intensity at large depth of 
matter, muon-induced hadron flux underground, background produced by 
atmospheric muons in underground neutrino experiments, etc. 

Several models [1-4] have been developed to describe the muon-nuclear in- 
elastic interactions, nonetheless uncertainties of this process are much larger 
than for purely electromagnetic interactions. The reason is that the bulk of 
this process is characterized by low squared four- momentum transfer Q 2 . The 
smallness of the Q 2 does not allow to use of the perturbative QCD (pQCD) 
for calculation of nuclear structure function (SF) and phenomenological mod- 
els such as the Regge or General Vector Dominance Model (GVDM) have to 
be used. The parameterization of the nucleon SF obtained in these models, 
depends on free parameters which can be determined from a fit of experi- 
mental data and can be applied in the range of Q 2 < 1 — 3 GeV 2 . This range 
is often referred to as photo-production. However, they fail to describe deep 
inelastic scattering (DIS) data at high Q 2 . The pQCD (NLO QCD) gives a 
good description of the structure functions at Q 2 > 3 GeV 2 . So, a model 
to combine various aspects of these approaches is needed to describe the Q 2 
behavior of nucleon SFs over complete range from photo-production to DIS. 

The widely used approximation [4] of the muon photo-nuclear cross- 
section was obtained twenty years ago in the framework of GVDM. Exper- 
imental data in the region of Q 2 < 100 GeV 2 and x > 0.01 were used for 
determing the parameters. Recently, precise data [5] on SF in wide ranges 
of Q 2 (0.06 < Q 2 < 10 4 GeV 2 ) and x (10~ 6 < x < 0.98) have been obtained 
and new nuclear effects (anti-shadowing and EMC-effect) were observed. 

The main goal of this work is the calculation of the muon-nucleus inelastic 
cross-section, based on the modern nucleon SF and on the present knowledge 
of nuclear effects. 

The paper is organized as follows. In Sec. 2 we give the general relations 
and definitions used for description of neutral current charged lepton-nucleon 
scattering. The procedure of the calculation of nucleon SF, using CKMT 
Regge model CKMT [6] and MRS99 parton distribution function (PDF) [7] 



2 



is described in Sec. 3. In Sec. 4, nuclear effects and their parameterization are 
described. The total cross-section, the muon energy losses, and the angular 
distributions of scattered muons are given in Sec. 5. In the conclusion we 
summarize the main results of the paper. 

2 Neutral current charged lepton-nucleon 
scattering cross-section 

The cross-section of neutral current charged lepton scattering off nucleon 

l(k) + N(p)^l(k')+X(p') (1) 

is given by the sum of contributions from the processes as shown in Fig.l Here 
k(E, k) and k'(E', k') are initial and final lepton four-momenta, q = k — k' is 
the virtual photon and Z-boson momentum, p and p' are the initial nucleon 
momentum and the total momentum of the final hadrons X, respectively. 
This process can be described by the four-momentum transfered Q 2 , Bjorken 
x and lepton's energy loss v (or inelasticity y) defined as 

Q 2 = - q 2 = (k - k') 2 , x = g, 
M p ■ k 



In laboratory system 



Q 2 = 2(EE' - W) -2m 2 , x Q 



2Mv 



v = E-E\ y=^ (3) 

where M and m are nucleon and lepton masses, respectively. 

The general form of the differential cross-section for scattering of charge 
non-polarized lepton on non-polarized nucleon, summed over the final lepton 
polarizations can be expressed as: 

d2(jl ~' l+ W f E^ + (x,Q 2 )-I^ + (x,Q 2 ) + 



dvdQ 2 Q 4 E 2 



x,Q 2 ) , (4) 
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where 



E i~,i + = 2x F?(x, g 2 )Fi + F?(x, Q 2 )Y 2 , (5) 

I l ~> l+ = P Z { gv(2xF((x,Q 2 )Y 1 + Fi(x,Q 2 )Y 2 ) T 

g A xFi(x,Q 2 )Y 3 y (6) 



^ + = P 2 z \{g 2 v + g 2 A )(2xF 1 z (x,Q 2 )Y 1 + F 2 z (x,Q 2 )Y 2 ) 

T2g v gAxF 3 z (x,Q 2 )Y 3 y (7) 



and 



Y 1 = (Q 2 - 2m 2 ) ^, F 2 = [2E(E - v) - 



91 il 

2 'Ji/ 



F 3 = (^-^). (8) 
Here the term P z accounts for the Z° propagator is 

z V2 2na Q 2 + m 2 x ' { ' 

where Gj \/2 is the Fermi constant, a = 1/137 is the fine structure constant, 
and m z is Z-boson mass. The constants of lepton weak coupling g v and g A 
are 

g v = -- + 2sm 2 9 w g A = --, (10) 

where 9w is Weinberg angle. 

The functions Ff' z are the electromagnetic (7-exchange) and neutral 
current (Z-exchange) structure functions, respectively. The functions F/ 
correspond to decomposition over invariant functions of the tensor (j-Z in- 
terference) 

W 1 ~J2{(p'\ J a\p)(p\ J p\P') + 

p'\jZ\p)(p\J$\p')}5(p'-p-q). (11) 
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The upper sign in Eq.(6) and Eq.(7) corresponds to lepton scattering 
(e~,/i~), and the lower sign is for antilepton(e + , scattering. 

The term proportional to the function F/, is due to interference between 
the electromagnetic scattering amplitude and the axial-vector current weak 
interaction amplitude. The amplitudes have opposite C parities, so the cor- 
responding terms have opposite sings for lepton and antilepton scattering. 
At low Q 2 the 7 — Z interference term is much smaller then the 7-exchange 
one, but it increases linearly with Q 2 (Eq.(9)) and becomes comparable with 
the 7-exchange term at Q 2 £sl0 3 GeV 2 . 

In terms of parton distributions in LO-approximation, the SFs can be 
written as 

F ? = \T,<Z (/, + /*)> (12) 

*?=*E e ; (/,+/*), (13) 
f i z =^eK 2 +°3 (/*+/*)> (14) 
f?=*t,{v1 (1 5 ) 

Fi = 2j2v q a q (f q -f,), (16) 

and 

^ = EV, (fq + fq), (17) 

F 2 ^2^e^(/ g + / 5 -), (18) 

Fi = 2'£e q a q (f q -f q ). (19) 

Here / 9 and f q are parton distribution functions in the proton, e q , v q and 
a q are the charge, vector and axial-vector week couplings of quarks. For 
up-quarks (u, c, t) they are 

2 1 4 - 1 

e u , c ,t = - u„ )C ,t = 2 - 3 sm *w au,c,t = g ( 20 ) 

and for down-quarks (d,s,b) 

1 1 2 2. 1 

ed, s ,fe = -- = -- + - sm 9 W a d:S:b = -- (21) 
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It is seen from Eqs.(4) and (9) that the main contribution to the total cross- 
section is due to photo-production (low Q 2 process) , however at fixed outgo- 
ing muon energy the large scattering angle corresponds to high Q 2 

cos6 = (EE' -Q 2 /2-m 2 ) / | k || k' | . (22) 

So, for calculation of muon scattering at large angles it is necessary to know 
the behavior of the nucleon SFs in the wide range of Q 2 ~ 0.01 — 10 6 GeV 2 . 

3 Low and high-Q 2 approximation of the nu- 
cleon structure functions 

At high-Q 2 the QCD predictions for the nucleon SFs are obtained by solving 
the DGLAP evolution at NLO approximation in MS or DIS schemes. These 
equations yield the parton distribution functions at all values of Q 2 provided 
that the PDF are given as a function of x at some input scale Qq=1.2-5 
GeV 2 . The latest global fits performed by several group (MRS99 [7], GRV98 
[8], and CTEQ5 [9]) give a good description of the experimental data. At Q 2 
below Ql perturbative QCD fails to describe data and phenomenological non- 
perturbative (GVDM or Regge) models are required. A considerable number 
of non-perturbative models have been developed [10]- [12] recently. These 
models predict a correct limit of F 2 at Q 2 =0 and give a good description of the 
SFs at low and medium Q 2 . Thus, neither the non-perturbative approaches 
nor pQCD can be expected to describe the Q 2 behavior of the SFs over 
the complete the range from photo-production to DIS. A number of models 
combining QCD and phenomenological approaches have been developed to 
describe data in the transition region of Q 2 (see review [13]). In this paper 
we use the CKMT model [6] at low and moderate Q 2 and the MRS99 fit of 
PDF [7] at high-Q 2 . 

The CKMT model proposes the following parameterization of the proton 
structure function 

Ff (*, Q 2 ) = F*{x, Q 2 ) + F p NS (x, Q 2 ) (23) 

The singlet term 

F$(x, Q 2 ) = A s x~^ (1 - x) n ^ \J^-j (24) 
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corresponds to the Pomeron contribution which determines the small-x be- 
havior of sea quarks and gluons. The dependence of effective intercept of the 
Pomeron, A on Q 2 is parametrized as 



A(Q 2 ) = A 1 + 



2Q 2 



Q 2 + dJ ' 

The x — > 1 behavior of F s (x, Q 2 ) is determined by the function 



(25) 



n(Q 



Q 2 + c, 



(26) 



The parameterization for the non-singlet term which corresponds to the sec- 
ondary (f,A 2 ) reggion (valent quark) contribution, is 

FUx, Q 2 ) = Bx^ (1 - xf^ l^f-^ ** , (27) 

where behavior at x — > is determined by the secondary reggion intercept 
ocr. The valence quark distribution can be separated into the contributions 
of the u and d valence quarks by replacing 



where 



F p NS {x, Q 2 ) = xU v (x, Q 2 ) + xD v (x, Q 2 ), 

-)2 \ a R 



xU v (x, Q 2 ) = B u x {l ' aR) (1 - x) n{Q2) ^ 



Q 



xD v (x, Q 2 ) = B d x (1 - aR) (1-x) 



n(Q 2 )+l 



Q 2 + bJ 

' Q 2 \ aR 



\Q 2 + b / 

Normalization conditions for valence quarks in proton 



(28) 

(29) 
(30) 



1 r 



X 



xU v (x,Q 2 ) 



dx = 2e. 



dx = e 2 d 



- [xD v (x,Q 2 ) 
fix the values of parameters B u and B d at Q 2 = Q\. 



(31) 
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The limit of Q 2 = corresponds to interaction of real photons. The total 
cross-section for real photons can be written as 



F2(x,Q 2 ) 



Q 2 =0 



One can see from Eqs.(23),(27), and (28) that F 2 « Q 2 at Q 2 
v. Thus, the parameterization 



(32) 
and fixed 



1_A °(2Mi/) A ° + 
(B tt + £ d )&- a " i (2Mi/) a *- 1 



(33) 



takes place in the CKMT model. 

In this way we find parameterizations of both F% and 'jp cross-section with 
7 free parameters: a, b, c, d, A , or, and Ag. To determine the parameters 
we have made a joint fit of the cflot 

data and NMC, E665, SLAC, ZEUS, 
and HI data on the proton SF F 2 in the region 0.11 < Q 2 < 5.5 GeV 2 
and 10~ 6 < x < 0.98 [5]. As initial condition for the values of different 
parameters, we used those obtained in the previous fit in Ref.[6]. A global fit 
results in the following values of parameters (all dimensional parameters are 
in GeV 2 ): a=0.2513, 6=0.6186, c=3.0292, d=1.4817, A =0.0988, a R =0.4056, 
and As— 0.12. The values of the parameters B u = 1.2437 and B d = 0.1853, 
were determined from normalization conditions for valence quarks (at Ql=2 
GeV 2 ). The quality of the description of all experimental data is very good 
and x 2 /d-o.f. =754. 8/600, where only the statistical errors have been used. 
Recently, a modified version of the CKMT model used new data on F$ at 
low Q 2 has been published [14]. The values of the main parameters are in a 
good agreement with those, obtained in the present work. 

So, for calculation of F 2 in the entire region of Q 2 we use the CKMT 
model at Q 2 < 5 GeV 2 , the MRS99 PDF at Q 2 > 6 GeV 2 and a linear fit 
between Ff (CKMT) and Ff (MRS99) in the region 5 < Q 2 < 6 GeV 2 . The 
result of the fit of F$ and <t*°' is shown in figures 2 (F$ vs. x for different value 
of Q 2 ) and 3 (F$ vs. Q 2 for different value of x) along with the experimental 
data [5] . The cross-section a^t as a function of W 2 = M 2 + 2Mv — Q 2 is 
shown in Fig. 4 (the data from Refs.[15],[16]). 

A good description of experimental data is obtained at all x and Q 2 
values. It should be noted that: 
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a) the recent ZEUS BPT97 data [17] were not included in our fit, but are 
in agreement with the CKMT model prediction at Q 2 < 0.1 GeV 2 , 

b) the rise of Ff at low x and low Q 2 is well described by the CKMT 
model with the slope A =0.0988 while the experimental value is 0.102±0.07 
[17], 

c) the <t*°* values found by the ZEUS collaboration are the result of a 
phenomenologically motivated extrapolation. 

In figures 2-4 the SF Ff and a^t that have been obtained by Bezrukov 
and Bugaev [4] and used for calculation of the muon photo-nuclear cross- 
section are shown also. At x < 10~ 3 the SFs rise slower than the present 
data indicate. An the other hand in the region x > 0.01 and Q 2 > 5 GeV 2 
the SFs are overestimated. 

The CKMT parameterization gives separate contributions of valence 
quarks, sea quarks, and gluons. We used this peculiarity for parameteri- 
zation of the neutron SF F 2 , that can be extracted from the deuteron F 2 
and the proton F 2 data, using the relation 

F 2 d = \ \Fl(x) + F?(x)\ , (34) 

and the Gottfried sum rule 

i i l 

Sg = J (Ff -F?)-£ = -f (d v - u v )dx --J (d-u)dx. (35) 



In the case of SU(2) symmetric sea, PDF of d equals to—u and therefore 
50=1/3. However, the NMC collaboration [5] gives ,Sg=0.235±0.026 at Q 2 = 
4 GeV 2 , i.e. significantly below 1/3 and shows that F 2 (x) — F£{x) — > and 
F 2 ra (x) / ' F 2 {x) — > 1 at x — > 0. Taking into account these results the singlet 
term of the F 2 has to be modified. Because of the isotopical invariance of 
strong interaction the non-singlet term F^ s is 

F^ 5 (x, Q 2 ) = ^xU v (x, Q 2 ) + AxD v (x, Q 2 ) (36) 

where xUv{x,Q 2 ) and xD v {x,Q 2 ) are given by Eq.(29) and Eq.(30). The 
singlet term 

/ n 2 \1+A(Q 2 ) 

FS(x, Q 2 ) = A s x~^ (1 - x) n ^ J (37) 



9 



has an additional free parameter r. The value of this parameter obtained 
from fit of the F 2 data [5] in the region Q 2 < 5. GeV 2 (all other parameters 
were fixed by the fit of the and a%) is r = 1.8152. The quality of the 
description of data is good with value of \ 2 /d.o.f = 611.1/453 for the SF F 2 
data and x 2 /d-o.f = 452.8/380 for the Fg/Fg data, where only statistical 
errors have been used. 

For calculation of the SF F 2 n in the entire region of Q 2 we used the 
approximation of Eq.(37) and Eq.(38) at Q 2 < 5. GeV 2 , the MRS99 PDF 
at Q 2 > 6 GeV 2 and linearly fit between F 2 n (CKMT) and F 2 n (MRS) in the 
transition region 5 < Q 2 < 6 GeV 2 . The SF F 2 d (Fig.5) and F 2 / F 2 (Fig.6) 
as function of Q 2 for different values of x are shown with experimental data. 
Fig. 7 shows the Fg - F 2 n vs. x at Q 2 = 4 GeV 2 . The calculations are m 
agreement with the NMC data [18]. 

For calculation of the cross-section of lepton-nucleon scattering it is nec- 
essary to know the behavior of SF 1xF\ in a wide range of Q 2 and x also. 
This SF can be expressed using the longitudinal SF F L 

( 4M 2 x 2 \ 

F L =ll + ^-\F 2 -2xF 1 . (38) 

Then 

1 / AM 2 x 2 \ , . 

where 

Fl 

R= (l+AM 2 x 2 /Q 2 )F 2 -F L (40) 

Perturbative QCD describes reasonably well the available data on the 
ratio R(x, Q 2 ) at large values of Q 2 and very little is known about possible 
extrapolations towards the region of low Q 2 . In the limit of Q 2 — > the SF 
Fl has to vanish as Q 4 (for fixed v ) and therefore R ~ Q 2 . At x < 0.01 and 
Q 2 < 0.5 GeV 2 the experimental results are poor. Data show a small value 
of R at moderate values of x and a possible increase of R as x decreases. 
The data come from experiments carried out on different targets and the 
differences R A — RP are consistent with zero and do not exhibit any significant 
dependence on x [19]. 

Data can be fitted by a parameterization i?(SLAC98) [20]. However, this 
fit should not be used at Q 2 < 0.35 GeV 2 . In this region we used the GVDM 
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asymptotic of R at Q 2 — > given by [21] 

Rgvmd(Q 2 ,x) « 



Q 2 + m 2 ' 



where m p = 0.77 GeV is the p meson mass. Thus, at Q 2 > Ql — 1.4 GeV 2 
the function Q 2 ) is calculated as follows: 

R(r 2 ] _ j R(MRS99) at x < lO" 3 

J - | ^(^^gs) at x > 5.10-3 l 4i J 

In the region 10" 3 < x < 5 • 10" 3 a linear fit between i?(MRS99, x = 10~ 3 ) 
and ^(SLACgS, x = 5 • 10- 3 ) is used. At Q 2 < Q 2 , 

i?(a:,g 2 ) = i? G yMD(a:,g 2 ) = C(x) Q * 2 , (42) 

*V + m p 

where the function C(x) is determined by normalization condition at Qq 

Rgvmd(x,Q 2 )=R(x,Q 2 ), (43) 

and the function R(x,Q 2 ) ) is calculated using Eq.(40). Fig. 8 shows the ex- 
perimental values of R as a function of Q 2 in four ranges of x along with the 
result of parameterization Eq.(41)-Eq.(42). In the region of low Q 2 < Qq 
it decreases with Q 2 at all values x, but the dependence on x is not strong 
(Fig. 9). However, extrapolation of R outside the kinematical range of data: 
namely at Q 2 — > and x — > 0, based on the presently available data, is a 
rather delicate problem. 

In Figs. 10 and 11 we show the results of calculations of the differential 
cross-section of neutral current (NC) e ± p scattering da/dQ 2 and da/dy at 
high Q 2 . The cross-section da/dQ 2 decreases by six orders of magnitude 
between Q 2 = 400 GeV 2 and 4000 GeV 2 . This decrease is due to the photon 
propagator. The cross-section da/dy is shown for different Q 2 regions. For 
Q 2 > 400 GeV 2 the bulk of the cross-section is concentrated at small values of 
y. For Q 2 > 10 4 GeV 2 the differential cross-section is approximately constant 
with y. The predictions using the MRS99 PDF give a good description 
of measured cross-sections. NC scattering at high Q 2 is sensitive to the 
contribution due to the Z°. According to Eqs.(4)-(9) the Z° contribution 
reduces approximately by 25%(12%) for e + p (e~p) cross-section at Q 2 > 10 4 
GeV 2 . 
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4 Nuclear structure functions 



The SFs measured for different nuclei A are found to differ from the SF 
measured on deuteron [24,25]. The modifications are usually observed as 
a deviation from unity of the ratio r A l d = F A /F d , where F A and F d are 
the SFs per nucleon measured in a nucleus and in deuteron, respectively. 
Neglecting nuclear effects in the deuteron, F d can approximately stand for 
a isospin averaged nucleon SF, F^ = (F$ + F 2 n ) /2. Different nuclear effects 
are observed in different region of x. 

(i) Shadowing at x <0.1. The ratio r A l d is smaller than unity. The ex- 
perimental data cover the region x >10 -4 and r A l d decreases with decreasing 
x. Shadowing increases with nuclear mass A and weakly depends on Q 2 . 

(ii) Anti-shadowing at 0.1 < x < 0.2. The NMC data have established a 
small (a few percent) but statistically significant excess over unity. Within 
the accuracy of the data no significant Q 2 -dependence of this effect has been 
found. 

(iii) EMC effect at 0.2 < x < 0.8. The measured ratio r A/d decreases 
as x rises and has a minimum at x — 0.6. The magnitude of this depletion 
grows slowly with nuclear mass number. The data imply that a strong Q 2 - 
dependence of the r A l d is excluded in this region also. 

(iv) Fermi motion. At x > 0.8 the ratio r A ^ d rises above unity but exper- 
imental information is rather scarce. 

Investigations of differences between the longitudinal-to-transverse cross- 
section ratio R = olI ot for different nuclei showed that R Al — R Aj is com- 
patible with zero. This implies that nuclear effects influence both SFs F\ 
and i<2 in a similar way. 

At the moment, there is no unique theoretical description of these effects; 
it is believed that different mechanisms are responsible for them in different 
kinematical regions. For example, the EMC effect indicates that the averaged 
momentum carried by valence quarks in nuclei is reduced, relative to free 
nucleon. It has been shown in [26] that the pattern of the function r A / d (x) 
has a universal shape in the range 10~ 3 < x < 0.96 and in the range of nuclei 
mass A >4. Namely, the ratio F A (x)/ F d (x) could be well approximated with 
phenomenological functions in different regions of x. At x > 0.3 

r A / d (x) = l-m b (A)a osc (x), (44) 
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where the A dependence of m b can be approximated as 



m b (A) =M b (l- N S {A)/A) and M b = 0.437 (45) 

for A 7^4. The number of nucleons N S (A) at the nucleare surface is given by 
the Woods-Saxon potential 

1 + exp Ur — r (A))/a\ 

r (A) 

with values of parameters po=0.17 fm~ 3 , a=0.54 fm and 

r (A) = 1A2A 1/3 - 0.86A- 1/3 . (47) 

The function a osc (x) is 

«^*> = a -■^{(H)-'* (?-?)}■ (48) 

where u = 1 — x, c = 1 — x 2 , X2 =0.278, A=0.5 and \i = m n /M (m^ is pion 
mass). At 10~ 3 < x <0.3 the function is given by 

r A/d (x) =x mi (l + m 2 )(l-m 3 x), (49) 

with 

m i = M i (l-N s (A)/A), (50) 

where M^O.129, M 2 =0.456 and M 3 =0.553. We used Eq.(49) for calculation 
of r A/d up to x < 10 -3 . 

The value of function of A was obtained in the following way. 

The experimental data [26], show that in the region 5 • 10~ 3 < x < 0.1 the 
ratio r A l d decreases with x. Generally, small x corresponds to small Q 2 , 
therefore the approach of real photon interaction can be used. So, for x — > 0, 
r A l d — * i] A =a^A/ Aa^N, where <t 7 a is photon-nuclear cross-section and <t 7 at 
is the photon-nucleon cross section averaged over proton and neutron. The 
expression for the function r\ A has been obtained in Ref . [4] , using the optical 
nuclear model 

rj A = 0.75G(z) + 0.25, (51) 
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where the function G(z) is 



--l + e-(l + z) 



(52) 



and z = 0.00282A 1 / 3 <T 7A r(z/). Using Eq.(35) with the values of parameters 
obtained in this work, we can write the averaged photon-nucleon cross-section 
as follows 



a lN = ^{a 1P + a in ) = 112.2 (0.609z/ - 0988 + 1.037^ 5944 ) 

In the range i < 1, Eq.(49) reduces to 

r A ' d {x) =x mi (l + m 2 ). 
Then, from the asymptotic condition 

r A / d (x ) = 0.75G(z) + 0.25 
we obtain the expression for xq 



x = 



(0.75G(z) + 0.25) 



l/mi 



.1 + m 2 

At x < x we assumed that the function r A l d is constant and 

r A / d (x ). 



r A/d (x) 



(53) 

(54) 
(55) 

(56) 

(57) 



The results of approximation of the ratio r A l d are presented in Fig. 12 as a 
function of x for different nuclear targets and are in a good agreement with 
data. 

Taking into account the nuclear effects, the nuclear SFs Fi(x,Q 2 ) and 
total photon-nuclear cross section can be written as: 



F A (x,Q 2 )=Ar A / d (x,A)F»(x,Q 2 



(58) 



and 



a lA (u) = Aa lN (u)[0.75G(z) + 0.25]. (59) 

The calculated cross sections o" 7 yi(z/) are shown in Fig. 13 as functions of the 
real photon energy for nuclei C, Cu and Pb. 
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5 Muon inelastic scattering in standard rock 



We have calculated the main characteristics of inelastic muon scattering in 
standard rock ( A = 22, Z — 11 and p = 2.65 g/cm 3 ). The spectra of muon 
energy loss (Na v is Avogadro's number) in single interaction, 

Qmax 

N A v da^ A N Av ?" da 2 A 
A dv A J dvdQ 2 v ' 

Q 2 ■ 

are shown in Fig. 14 as functions of inelasticity y for different muon energies. 
The energy dependence of the total cross section, 

^a{E)= J ^dv, (61) 

and muon energy loss, 

- % / < 62 » 

are shown in Fig. 15. The allowed kinematical region for the variables v and 
Q 2 is determined by the following equations: 

Q 2 = 2 (EE'- | k 1 1 k' | cos 6) - 2m 2 (63) 

at cos# = ±1, and 

Q 2 = 2Mv + M 2 — W 2 . (64) 

The results from Ref.[4] are given in these figures for comparison also. It is 
necessary to note, that the cross-section and muon energy loss [4] have been 
calculated taking into account the shadowing effect only. The cross-section 
of inelastic muon scattering obtained in the present work is larger by factor 
1.2 and the muon energy losses b n (E) is larger too by ^8% at E =10 3 GeV 
and by ^30% at E =10 6 GeV. As a result, the total energy losses (the sum 
of bremsstrahlung, pair-production and inelastic muon scattering) increases 
by wl% at E =10 3 GeV and by ^4% at E =10 6 GeV. These differences are 
mainly due to contributions of small x and small Q 2 where the modern SFs 
are large than SFs used by Bezrukov and Bugaev. 
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The probabilities P(> 9, > v) of muon scattering in single interaction at 
the angles larger than 9 with outgoing muon energy E' > vE are shown in 
Fig. 16 as a function of 9 for different values of v and primary muon energies. 
The results are given for /i~ and ji + scattering. The main peculiarities of the 
inelastic muon scattering are: 

(i) At fixed values of 9 and E' the probability decreases very quickly 
with the initial muon energy E. For example, for 9 > 2° and E'=10 GeV, 
P=6.3-10" 4 at P=10 2 GeV and P =3.8-10~ 6 at P=10 3 GeV. 

(ii) At fixed values of 9 and E the probability increases with decreas- 
ing outgoing muon energy E '. For 9 >2° and P=10 3 GeV, P=7.8-1CT 7 at 
E' >10 2 GeV and P =3.8-l(T 6 at E' >10GeV. 

(iii) At fixed values of E and E' the mean values of x and Q 2 ((x) and 
(Q 2 )) increase with 9. For muon energies E=10 2 GeV and E' >0.1E the 
values of (x) and (Q 2 ) increase from (x)=0.12 (anti-shadowing region) and 
(Q 2 )=0.75 GeV 2 at 0=0.25° up to (x)=0.25 (EMC-region) and (Q 2 >=34 
GeV 2 at 9=6°. For energies E=10 3 GeV and E' >0.1E: (x)=0.09 (shad- 
owing region), < Q 2 >=28 GeV 2 at 0=0.25° and (x)=0A6 (anti-shadowing 
region) and (Q 2 )=925 GeV 2 at 0=6°. Thus, the probability of scattering at 
large angles is suppressed by the EMC-effect. 

6 Conclusions 

In this work we studied inelastic muon scattering off nuclei. 

1. It is shown that the combination of the SF F2, predicted by the CKMT 
model at low and moderate values of Q 2 , and the MRS99 PDF at high Q 2 
gives a good description of the data over the complete measured region from 
photo-production to DIS. In particular, the CKMT model well describes the 
rise of the F$ at low x and Q 2 with the slope Ao = 0.0988. Furthermore, 
in the framework of this model the expression for the neutron SF F£ can be 
obtained. The result is in a good agreement with the F 2 /F 2 and F 2 — F 2 
data. 

2. The MRS99 PDF well describe the differential cross-sections da/dQ 2 
and da/dy, calculated taking into account not only electromagnetic current 
contribution but also contributions of neutral current and 7 — Z interference. 
The 7 — Z interference contribution is clearly seen at high Q 2 > 10 3 GeV 2 . 

3. The nuclear effects modify the nucleon SFs in the entire measured 
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region of x and Q 2 . The modification depends very slowly on Q 2 and increases 
with A. 

4. The obtained SFs have been used for calculations of the muon-nucleus 
scattering cross-section, muon energy losses and muon angular distributions 
in inelastic interaction taking into account nuclear effects and contributions of 
neutral current and 7 — Z interference. As a result, the total cross-section and 
energy losses increase with muon energy faster then predicted in Ref . [4] . The 
scattering of high energy muons (E > 10 3 ) GeV at large angles is suppressed 
by EMC-effect and 7 — Z interference. 
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Fig. 2a: The proton structure function at low Q 2 as a function of x. 
Solid lines are our calculations using the CKMT model and dashed lines 
calculations of Bezrukov and Bugaev [4]. Points are results of experiments 
ZEUS and E665 [5,17]. 
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Fig. 2b: The same as Fig. 2a but at moderate Q 2 . 
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Fig. 3: The proton structure function Ff as a function of Q 2 at fixed values of 
x. The solid curves have been obtained by using the CKMT+MRS99 model 
and dashed curves are result from Ref.[4]. Data set is due to ZEUS, HI, 
E665, NMC, and SLAC experiments [5]. For clarity an amount Cj=13.6-0.6i 
is added to where i=l(16) for the lowest (highest) x value. 
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Fig. 4: Total cross-section a^t as a function of W 2 . Solid line is our calcula- 
tions using combined (CKMT+MRS99) model and dashed one is from work 
[4]. Experimental data are from [15] at low energies and from ZEUS and HI 
[16] at higher energies. 
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Fig. 5: The deutron structure function F% as a function of Q 2 at fixed values 
of x. Data set is due to BCDMS, E665, NMC, and SLAC experiments [5]. 
For clarity an amount Cj=13.6-0.6i is added to _F 2 d where 2=1(10) for the 
lowest (highest) x value. 
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Fig. 6: The ratio of structure functions F^jF^ as a function of Q 2 at fixed 
values of x, as compared to NMC data [5]. For clarity an amount Cj = 
13.6 — 0.6i is added to F%fF% where for the lowest (highest) x value. 
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Fig. 7: The difference Ff — F 2 n at Q 2 =4GeV 2 as a function of x, as compared 
to NMC data [5]. 
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Fig. 8: The ratio R(x, Q 2 ) as a function of Q 2 at fixed x. Data are from 
Refs.[5] and [20]. 
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Fig. 9: The ratio R(x,Q 2 ) as a function of x at fixed Q 2 . Data are from 
Refs.[5] and [20]. 
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Fig. 10: The differential cross-section da/dQ 2 of neutral current proton scat- 
tering off electron (positron) as a function of Q 2 . Points are experimental 
results of ZEUS [22] and HI [23] experiments. 
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Fig. 11: The differential cross-section da/dy of neutral current proton scat- 
tering off electron (positron) as a function of y. Experimental data are due 
to ZEUS [22] experiment. 
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Fig. 12: The comparison of ratio F 2 A /F 2 d , measured by SLAC, BCDMS and 
NMC groups and approximations given in work [26]. 
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Fig. 13: Total cross-section <t 7 a(^) for C, Cu, and Pb as a function of real 
photon energy. Results of our calculations (solid curves) and calculations by 
Bezrukov and Bugaeev [4] (dashed curves) as compered to experimental data 
[27]. 
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Fig. 14: The spectra of muon energy loss due to muon inelastic scattering 
in standard rock as a function of y for fixed muon energies, as compered to 
calculations by Bezrukov and Bugaev [4]. (y is the fraction of lost by 
muon in single interaction.) 
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Fig. 15: The total cross-section g^a and muon energy loss b n for muon inelastic 
scattering in standard rock as a function of muon energy E^, as compared 
to those from Ref.[4]. 
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Fig. 16: The probabilities P(> 6, > v) of muon scattering in single inter- 
action at angle larger than 6 with energy of outgoing muon E' > vE as a 
function of 9 for fixed values of v and incoming muon energies. Solid (dashed) 
lines are for positiv (negativ) muon charge. 
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